Fully Characterizing Axially Symmetric Szekeres 
Models With Three Data Sets 



December 17, 2012 



Marie-Noelle Celerier 1 a , Priti Mishra 2 b and Tejinder P. Singh 2 c 

1 Laboratoire Univers et THeories (LUTH), Observatoire de Paris, CNRS, Universite 
Paris-Diderot, 5 place Jules Janssen, 92190 Meudon, France 

2 Department of Astronomy and Astrophysics, Tata Institute of Fundamental Re- 
search, Homi Bhabha Road, Colaba, Mumbai, 400005, Maharashtra, India 

e-mail: a marie-noelle.celerier@obspm.fr 
b priti@tifr.res.in 
c tpsingh@tifr.res.in 



Abstract 

Inhomogeneous exact solutions of General Relativity with zero cosmo- 
logical constant have been used in the literature to challenge the ACDM 
model. From one patch Lemaitre-Tolman-Bondi (LTB) models to axially 
symmetric quasi-spherical Szekeres (QSS) Swiss-cheese models, some of 
them are able to reproduce to a good accuracy the cosmological data. It 
has been shown in the literature that a zero A LTB model with a central 
observer can be fully determined by two data sets. We demonstrate that 
an axially symmetric zero A QSS model with an observer located at the 
origin can be fully reconstructed from three data sets, number counts, 
luminosity distance and redshift drift. This is a first step towards a fu- 
ture demonstration involving five data sets and the most general Szekeres 
model. 



1 The Szekeres solution 

We are interested here in finding an algorithm for building an axially symmetric 
quasi-spherical Szekeres (QSS) model from a set of observables. The problem to 
be solved is to determine the number of observables needed and to give a recipe 
allowing us to fully define the model from them. This problem has already been 
solved for Lemaitre-Tolman-Bondi models with a central observer by Mustapha 
et al. pQ. 

The QSS line element in comoving and synchronous coordinates is 
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where $ is a function of the t and r coordinates and 



S 



E = 



The Einstein equations give 
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where a dot denotes differentiation with respect to t and a prime, differentiation 
with respect to r. Equation ([3]) can be integrated as 
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Here, S(r),P(r),Q(r),k(r),M(r),tB(r) are all functions of the r coordinate. 



2 Set of differential equations 



In the axially symmetric case, we can set P = Q = [2J. The model is thus 
defined by three independent functions of r, plus a coordinate choice. We use 
this freedom to define the r coordinate such that, on the past light cone of the 
observer located at (t = to, r = 0), 
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where a hat denotes quantities computed on the past light cone. Substituting 
our coordinate choice into the total derivative of $ and using Eq. ([3]) with A = 0, 
we obtain 
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Now, we substitute Eq. into Eq. (jU), using Eq. ©, which gives 
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Then, we consider the luminosity distance on the axially directed past null 
geodesic, Dl{z) = Da{z) (1 + z) 2 which we substitute into the equation for Da 
as given by Bolejko and Celerier [3] and obtain a differential equation of the 
form 

F U 2 Dl/Az 2 ,dD~ L /dz, D~ L ,d 2 $/dr 2 ,d$/dr^,dp/dz,p,z,M,k,k',k") = 0, 
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The third differential equation needed proceeds from the equation for the red- 
shift and is of the form 
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We have obtained a set of three coupled differential equations for the un- 
known functions M, ft, and $, in which z and -Dl,P and their derivatives with 
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respect to z are observables. Now, to be able to determine the model, we need 
M, k, and $ as functions of r. Since the data are given in terms of the redshift 
z, we have to transform them into functions of r. We use thus the redshift 
drift, z'o(z) = (6z/Sto)(z) as given by Mishra et al. [3] and obtain a differential 
equation for z(r): 



We solve numerically the above set of differential equations with the origin 
conditions M(r = 0) = k(r = 0) = $(r = 0) = z(r = 0) = 0, and obtain 
M{r),k(r) then (E'/E)(r) and thus S(r), from the relation E'/E = S'/S, valid 
on an axially directed geodesic. With P(r) = Q(r) = and our coordinate 
choice, this is sufficient to fully determine the model. 



To obtain the functions, M, k, E and ts from the galaxy number count, n(z), 
supernova luminosity distance and redshift drift observations, we propose the 
following algorithm: 

f . Make a rigorous processing of the data to obtain smooth analytic functions 
of Dl(z),zq(z) and p{z) = mn(z), assumed to obtain on the axial null 
geodesic directed towards the observer located at the origin. 

2. Determine z{r) by solving numerically Eq. (jf f p . 

3. Determine 3?(r), M(r) and k(r) from the set of three coupled differential 
equations (|5])-([TD]). 

4. Knowing M(r),k(r) and determine E(r) (equivalently S(r)) by in- 
tegrating Eq. ([7]). 

5. Knowing now the sign of E, ts(r) follows from the resolution of the para- 
metric equations issued from Eq. ([5]) [5]. 

4 Conclusion 

We have given here a set of equations and an algorithm able to fully determine 
an axially symmetric QSS model of universe from three data sets assumed to be 
measured on an axially directed null geodesic issued from an observer located 
at the origin. This algorithm can only be applied to an ideal universe, where 
the redshift is monotonically increasing and where no shell-crossing is present 
to close the past null cone. However, since axial symmetry is not a very credible 
feature of our observed Universe, we do not claim this work to be more than a 
mere exercise designed to go beyond spherical symmetry and which will have to 
be generalized in the future to more physical Szekeres models. 
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